Generalizations of adjoint networks techniques for RLC interconnects model-order reductions

ABSTRACT

The adjoint network reduction technique has been shown to reduce 50% of the computational complexity of constructing the congruence transformation matrix. The method was suitable for analyzing the special multi-port driving-point impedance of RLC interconnect circuits. This paper extends this technique for the general circumstances of RLC interconnects. Comparative studies among the conventional methods and the proposed methods are also investigated. Experimental results will demonstrate the accuracy and the efficiency of the proposal method.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to generalizations of adjoint networks techniques, and more particularly to generalizations of adjoint networks techniques for RLC interconnects model-order reductions.

2. Description of Related Art

With the advances in deep submicron semiconductor techniques, the parasitic effects of interconnects has no longer been ignored [reference 2]. To completely characterize the signal integrity issues, interconnects are often modeled as large RLC networks. To this end, in order to reduce the computational time of the large-scale RLC interconnect networks, model-order reduction methods have been emerged recently [references 3, 6 and 9].

A consensus has emerged that of many model-order reduction techniques, the moment matching approach, or the so-called Krylov subspace projection method, seems to be the most viable one [references 3, 4, 5, 8 and 11]. In general, these methods can be divided into two categories: one-sided projection methods [references 8 and 11] and twosided projection methods [references 4 and 5]. The one-sided projection methods use the congruence transformation to generate passive reducedorder models while the two-sided ones can not be guaranteed. In recent works, we have proposed the adjoint network reduction technique to further reduce the cost about yielding the congruence transformation matrix [reference 7]. The method was suitable for analyzing the special multi-port driving-point impedance of RLC interconnect circuits.

The purpose of this paper is to extend the adjoint network technique for general RLC interconnect networks. First, relationships between an original MNA network and its corresponding adjoint MNA network are explored. Second, the congruence transformation matrix can be constructed by using the resultant biorthogonal bases from the Lanczos-type algorithms. Therefore, less storage and computational complexity are required in our proposed technique.

SUMMARY OF THE INVENTION

The main objective of the present invention is to provide improved generalizations of adjoint networks techniques for RLC interconnects model-order reductions.

The adjoint network reduction technique has been shown to reduce 50% of the computational complexity of constructing the congruence transformation matrix. The method was suitable for analyzing the special multi-port driving-point impedance of RLC interconnect circuits. This paper extends this technique for the general circumstances of RLC interconnects. Comparative studies among the conventional methods and the proposed methods are also investigated. Experimental results will demonstrate the accuracy and the efficiency of the proposal method.

Further benefits and advantages of the present invention will become apparent after a careful reading of the detailed description with appropriate reference to the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a mech circuit diagram;

FIG. 2 shows the frequency responses of the voltage Vout in FIG. 1; and

FIG. 3 shows the relative errors in the frequency responses in FIG. 2.

DETAILED DESCRIPTION OF THE INVENTION

The dynamics of RLC interconnect networks can be represented by the following modified nodal analysis (MNA) formulae [reference 8 and 11]: $\begin{matrix} {{{{M\frac{\mathbb{d}{x(t)}}{\mathbb{d}t}} = {{- {{Nx}(t)}} + {{Bu}(t)}}},{{y(t)} = {D^{T}{x(t)}}}}{where}{{M = \begin{bmatrix} C & 0 \\ 0 & L \end{bmatrix}},{N = \begin{bmatrix} G & E \\ {- E^{T}} & R \end{bmatrix}},{{x(t)} = \begin{bmatrix} {v(t)} \\ {i(t)} \end{bmatrix}},}} & (1) \end{matrix}$ x(t)εR^(n) is the state vector, u(t)εR^(m) is the input vector, y(t)εR^(p) is the output vector, and M,NεR^(nxn), BεR^(nxm),and DεR^(nxp) are so-called the MNA matrices. M and N containing capacitances in C, inductances in L, conductances in G and resistances in R are positive definite, M is symmetric and N is non-symmetric. E presents the incident matrix for satifying Kirchhoff's current law. x(t) contains node voltages V(t)εR^(nv) and branch currents of inductors i(t)εR^(ni), where n=n_(v)+n_(i). If the m-port driving-point impedance is concerned, then p=m and D=B.

Let the signature matrix be defined as S=diag(I_(nv),−I_(ni)) so that the symmetric properties of the MNA matrices are as follows [reference 11]: S ⁻¹ =S,SMS=M, and SNS=N ^(T)   (2) Under this situation, if port driving-point impedance is concerned, that is, each port is connected with a current source, then B^(T)=└B_(v) ^(T)0┘, where B_(v)εR^(nvxm), and SB=B.

The transfer functions of the state variables and of the outputs are defined as X(s)=(n+sM)⁻¹ and Y(s)=D^(T)X(s). Given a frequency s₀εC, let A=−(N+s₀M)⁻¹M and R=(N+s₀M)⁻¹B, where N+s₀M is assumed nonsingular. The Taylor series expansion of X(s) about s₀ is given by ${{X(s)} = {\sum\limits_{i = 0}^{\infty}{{X^{(i)}\left( s_{0} \right)}\left( {s - s_{0}} \right)}}},{where}$ X^((i))(s₀) = A^(i)R is the ith-order system moment about s₀. Similarly, the ith-order output moment about s₀ is calculated asY^((i))(s₀)=D^(T)X^((i))(s₀).

Suppose that the above system is large-scale and sparse. The model-order reduction problem is to seek a q-order system, where q<<n, such that $\begin{matrix} {{{\hat{M}\frac{\mathbb{d}{\hat{x}(t)}}{\mathbb{d}t}} = {{{- \hat{N}}{\hat{x}(t)}} + {\hat{B}{u(t)}}}},{{\hat{y}(t)} = {{\hat{D}}^{T}{\hat{x}(t)}}}} & (3) \end{matrix}$ where {circumflex over (x)}(t)εR^(q),{circumflex over (M)},{circumflex over (N)}εR^(qxm),{circumflex over (B)}εR^(qxm), and {circumflex over (D)}εR^(qxp). The corresponding ith-order system moment and output moment about s₀ is {circumflex over (X)}^((i))(s₀)=(−({circumflex over (N)}+s₀{circumflex over (M)})⁻¹{circumflex over (M)})^(i)({circumflex over (N)}+s₀{circumflex over (M)})⁻¹{circumflex over (B)} and Ŷ^((i))(s₀)={circumflex over (D)}^(T){circumflex over (X)}^((i))(s₀). The purpose of the moment matching is to construct a reduced-order system such that Ŷ^((i))(s₀)=Y^((i))(s₀) for 0≦i≦k−1, where k is the order of moment matching.

One conventional solution for moment matching is using the one-sided Krylov subspace projection method [references 3, 8 and 11]. The kth-order block Krylov subspace is defined as K(A,R,k)=colsp([RAR . . . A^(k−1)R]). K(A,R,k) is indeed the subspace spanned by X^((i))(s₀) for 0≦i≦k−1. The projection can be achieved by constructing V_(q)εR^(nxq),q≦km, from the Krylov subspace K(A,R,k). Under this framework, we have x(t)=V_(q){circumflex over (x)}_(q)(t) and the reduced-order model can be expressed as {circumflex over (M)}=V _(q) ^(T) MV _(q) ,{circumflex over (N)}=V _(q) ^(T) NV _(q) ,{circumflex over (B)}=B _(q) ^(T) B,{circumflex over (D)}=V _(q) ^(T) D   (4) It has been shown that X^((i))(s₀)=V_(q){circumflex over (X)}^((i))(s₀) and Ŷ^((i))(s₀)=Y^((i))(s₀) for 0≦i≦k−1. The reduced-order model is guaranteed stable {circumflex over (M)} and {circumflex over (N)} are positive definite. Furthermore, it will be passive if the multi-port driving-point impedance is concerned.

Two types of algorithms can be employed to generate V_(q) from the Krylov subspace: the Arnoldi type [references 8 and 11] and the Lanczos type [references 4 and 5]. We use the notation V_(q(A)) to denote the orthonormal basis generated from the block Arnoldi algorithm from the Krylov subspace K(A,R,k). Similarly, we use the notation V_(q(L)) and W_(q(L)) to represent the biorthogonal bases yielded from the block Lanczos algorithm from the Krylov subspaces K(A,R,k) and K(A^(T),D,k), respectively. In this case, W_(q(L)) ^(T)V_(q(L))=Δ_(q), where Δ_(q) is a full rank diagonal matrix. In the past, either V_(q(A)) or V_(q(L)) has been used to generate the reduced-order model (4).

Traditionally, W_(q(L)) and V_(q(L)) are used to perform oblique projection {circumflex over (M)}=−W _(q(L)) ^(T) AV _(q(L)) ,{circumflex over (N)}=W _(q(L)) ^(T)(I+s ₀ A)V _(q(L)) ,{circumflex over (B)}=W _(qT(L)) ^(T) R,{circumflex over (D)}=V _(q(L)) ^(T) D Although the reduced-order system can match up to 2k-order moments, this model can not be guaranteed to be stable and passive. Variations of the Lanczos-type algorithms have also been proposed. For example, if the m-port driving-point impedance is concerned (D=B), the symmetric block Lanczos algorithm has been investigated. In this case, W_(q)=D_(q)V_(q(L)), where D_(q) is a diagonal matrix, only a half of the computational cost and storage are required [references 3 and 5].

Using both W_(q(L)) and V_(q(L)) in the one-sided projection method is still possible. It can be achieved by the adjoint network reduction technique. The details can be developed in the following section.

One technique for model-order reductions is tot apply the corresponding adjoint MNA formulae. The adjoint network (or the dual system [reference 11]) of the system (1) is represented as [reference 10] $\begin{matrix} {{{M\frac{\mathbb{d}{x_{a}(t)}}{\mathbb{d}t}} = {{{- N^{T}}{x_{a}(t)}} + {{Du}(t)}}},{{y_{a}(t)} = {B^{T}{x_{a}(t)}}}} & (5) \end{matrix}$ The system transfer function and its ith-order system moment about s₀ are defined as X_(a)(s)=(N^(T)+sM)⁻¹D and X_(a) ^((i))(s₀), respectively. By using the information of X_(a) ^((i))(s₀), the following theorem has been shown in [reference 11]. Theorem 1: If a matrix U is chosen as the congruence transformation matrix such that X ^((i))(s ₀),X _(a) ^((j))(s ₀)}εcolsp(U) for 0≦i≦k−1 and 0≦j≦l −1,   (6) then Ŷ^((i))(s₀)=Y^((i))(s₀) for 0≦i≦k+l −1.

The technique can overcome the numerical instability problem when generating the basis matrix U if order k+l is extremely high. In this section, we summarize some properties of the adjoint network reduction technique. Relationships between X_(a) ^((i))(s₀) and X^((i))(s₀) are derived explicitly. It can be contributed to reduce the computational cost of constructing U.

In this subsection, we will show that the congruence transformation matrix U can be constructed with the resultant biorthogonal bases V_(q(L)) and W_(q(L)) from the Lanczos-type algorithms. The theorems will be disclosed as below.

By changing the state variables: x_(a)(t)=(N^(T)+s₀M)⁻¹z_(a)(t). Equation (5) can be rewritten as ${{{M\left( {N^{T} + {s_{0}M}} \right)}^{- 1}\frac{\mathbb{d}{z_{a}(t)}}{\mathbb{d}t}} = {{{- {N^{T}\left( {N^{T} + {s_{0}M}} \right)}^{- 1}}{z_{a}(t)}} + {Du}}},{{y_{a}(t)} = {{B^{T}\left( {N^{T} + {s_{0}M}} \right)}^{- 1}{{z_{a}(t)}.}}}$ The corresponding system transfer function is Z _(a)(s)=(N ^(T) +s ₀ M)(N ^(T) +sM)⁻¹ D=(N ^(T) +s ₀ M)X _(a)(s) . Thus the ith-order system moments of Z_(a)(s) and X_(a)(s) about s₀ can be derived as follows: Z _(a) ^((i))(s ₀)=(N ^(T) +S ₀ M)X _(a) ^((i))(s ₀)   (7) Through the introduction of Z_(a) ^((i))(s₀), the relationships between X_(a) ^((i))(s₀) and W_(q(L)) are observed in the following lemma. Lemma 1: Suppose that K(A^(T),D,k)=colsp(W_(q(L))). We have Z_(a) ^((i))(s₀)εcolsp(W_(q(L))),   (8) X_(a) ^((i))(s₀)εcolsp((N^(T)+S₀M)⁻¹W_(q(L)),   (9) for Osisk-1. Proof: First, from (7), we have $\begin{matrix} {{Z_{a}^{(i)}\left( s_{0} \right)} = {{\left( {N^{T} + {s_{0}M}} \right)\left\lbrack {{- \left( {N^{T} + {s_{0}M}} \right)^{- 1}}M} \right\rbrack}^{i}\left( {N^{T} + {s_{0}M}} \right)D}} \\ {= {\left\lbrack {- {M\left( {N^{T} + {s_{0}M}} \right)}^{- 1}} \right\rbrack^{i}D}} \\ {= {\left( A^{T} \right)^{i}{D.}}} \end{matrix}$ Thus colsp([Z_(a)⁽⁰⁾(s₀)Z_(a)⁽¹⁾(s₀)  …  Z_(a)^((k − 1))(s₀)]) = K(A^(T), D, k) = colsp(W_(q(L))). Second, the above result implies that Z_(a) ^((i))(S₀)=W_(q(L)){circumflex over (Z)}_(a) ^((i))(s₀) for 0≦i≦k−1. Thus (9) can be proven as below: X _(a) ^((i))(s ₀)=(N ^(T) +s ₀ M)⁻¹ W _(q(L)) {circumflex over (Z)} _(a) ^((i))(S ₀), for 0≦i≦k−1. Each X_(a) ^((i))(s₀) exists in the subspace spanned by columns of (N^(T)+S₀M)⁻¹W_(q(L)). This completes the proof. Theorem 2: Suppose that X^((i))(s₀)εcolsp(V_(q(L))), for 0≦i≦k−1, is a set of moments of x(s) about s₀. X_(a) ^((i))(s₀)εcolsp((N^(T)+s₀M)⁻¹W_(q(L))) for 0≦i≦k−1. Suppose that V_(q(L)) and W_(q(L)) are biorthogonal matrices generated by the block Lanczos algorithm. Let U=[V_(q)(N^(T)+s₀M)⁻¹W_(q)] be the congruence transformation matrix for model-order reductions in (4), then Ŷ(s ₀)=Y ^((i))(s ₀), for 0≦i≦2k−1,   (10) Proof: From the projection theory in Section 2 and Lemma 1, we have X ^((i))(s ₀)=U{circumflex over (X)} ^((i))(s ₀), for 0i≦k−1, and X _(a) ^((i))(s ₀)=U{circumflex over (X)} _(a) ^((i))(s ₀), for 0≦i≦k−1. Then the result can be shown by Theorem 1.

Theorem 2 demonstrates that a stable reduced-order model for general RLC circuits can be generated by the Lanczos-type algorithms. Moreover, 2k-order output moments are matched by performing k iterations of the algorithm.

Theorem 3: Suppose that X^((i))(s₀)εcolsp(V_(q(A))) for 0≦i—k−1 is a set of moments of X_((s)) about s₀. V_(q(A)) is theorthonormal matrix generated by the block Arnoldi algorithm. Let U=[V_(q(A))SV_(q(A))] be the congruence transformation matrix for model-order reductions in (4), then Ŷ ^((i))(s ₀)=Y ^((i))(s ₀), for 0≦i≦2k−1   (11) Proof: Let P be a q×n diagonal matrix. Since W_(q(L))=P(N^(T)+s₀M)V_(q(L)) and V_(q(L))=V_(q(A)) [references 3 and 5], it can be shown that the subspace spanned by U=[V_(q(A))SV_(q(A))] and U=[V_(q(L))(N^(T)+s₀M)⁻¹W_(q(L))] are the same.

A mesh twelve-line circuit, presented in FIG. 1, is studied to show the efficiency of the proposed method. The line parameters are resistance: 1.0 O/cm, capacitance: 5.0 pF/cm, inductance: 1.5 nH/cm, driver resistance 3 O, and load capacitance: 1.0 pF. Each line is divided into 50 sections. Therefore, the dimension of the MNA matrices is n=1198, m=1, and p=1, and D≠B. We set the expansion frequency s₀=1 GHz, the iteration number k=15, and use a total of 1001 frequency points distributed uniformly between the frequency range {0,15 GHz} for simulations. The frequency responses of the original model and the reduced-order model generated by the following projections: (1) W_(q(L)) and V_(q(L)); (2) V_(q(A)); (3) U=V_(2q(A)); (4) U=[V_(q(L))(N^(T)+s₀M)⁻¹W_(q(L)))]I are compared in FIG. 2. Their corresponding relative error, |Y(s)−Ŷ(s)|/|Y(s)|, are illustrated in FIG. 3.

The program was implemented in Matlab 6.1 with Pentium IV 2.8 GHz CPU and 1024 MB DRAM. The time to generate each reduced-order model and the corresponding average 1-norm relative error are summarized. The average 1-norm relative error is approximated by $\left( {\sum\limits_{i = 1}^{1001}{{{{Y\left( s_{i} \right)} - {\hat{Y}\left( s_{i} \right)}}}/{{Y\left( s_{i} \right)}}}} \right)/1001.$ Note that only 60% work is needed to generate the similar frequency response by using the proposed method.

Although the invention has been explained in relation to its preferred embodiment, it is to be understood that many other possible modifications and variations can be made without departing from the spirit and scope of the invention as hereinafter claimed.

REFERENCE

[1] Z. Bai, J. Demmel, J. Dongarra, A. Ruhe, and H. van der Vorst (edc.), Templates for the Solution of Algebraic Eigenvalue Problems: A Practical Guide, SIAM, Philadelphia, Pa., 2000.

[2] M. Celik, L. T. Pileggi, and A. Odabasioglu, IC Interconnect Analysis, Kluwer Academic Publisher, 2002.

[3] R. W. Freund, “Krylov-subspace methods for reduced-order modeling in circuit simulation,” Journal of Computational and Applied Mathematics, vol. 123, pp. 395-421, 2000.

[4] P. Feldmann and R. W. Freund, “Efficient linear circuit analysis by Pad_approximation via the Lanczos process,” IEEE Trans. On Computer-Aided Design of Integrated Circuits and Systems, vol. 14, no. 5, pp. 639-649, 1995.

[5] P. Feldmann and R. W. Freund, “The SyMPVL algorithm and its applications to interconnect simulation,” Proc. 1997 International Conference on Simulation of Semiconductor Processes and Devices, pp. 1113-116, 1997.

-   [6] S. Y. Kim, N. Gopal, and L. T. Pillage, “Time-domain macromodels     for VLSI interconnect analysis,” IEEE Trans. on Computer-Aided     Design of Integrated Circuits and Systems, vol. 13, no. 10, pp.     1257-1270, 1994.

[7] H. J. Lee, C. C. Chu, and W. S. Feng, “Interconnect modeling and sensitivity analysis using adjoint networks reduction technique,” in Proc. the 2003 International Symposium on Circuits and Systems, pp. 648-651, Bangkok, Thailand, 2003.

[8] A. Odabasioglu, M. Celik, and L. T. Pileggi, “PRIMA: passive reduced-order interconnect macromodeling algorithm,” IEEE Trans. on Computer-Aided Design of Integrated Circuits and Systems, vol. 17, no. 8, pp. 645-653, 1998.

[9] L. T. Pillage and R. A. Rohrer, “Asymptotic waveform evaluation for timing analysis,” IEEE Trans. on Computer-Aided Design of Integrated Circuits and Systems, vol. 9, no. 4, pp. 352-366, 1990.

[10] J. Vlach and K. Singhal, Computer Methods for Circuit Analysis and Design, 2nd ed., Van Nostrand Reinhold, 1993.

[11] J. M. Wang, C. C. Chu, Q. Yu, and E. S. Kuh, “On projection based algorithms for model order reduction of interconnects,” IEEE Trans. on Circuits and Systems-I: Fundamental Theory and Applications, vol. 49, no. 11, pp. 1563-1585, 2002 

1. Generalizations of adjoint networks techniques for RLC interconnects model-order reductions, the congruence transformation matrix U for the one-sided projection being defined as U=[V _(q,L) (N ^(T) +s ₀ M)⁻¹ W _(q,L)], wherein A. N and M matrixes are the matrixes of modified nodal analysis network of the system, such that ${{M\frac{\mathbb{d}{x(t)}}{\mathbb{d}t}} = {{- {{Nx}(t)}} + {{Bu}(t)}}},{{y(t)} = {D^{T}{x(t)}}}$ where x(t)εR^(n) is the state vector, u(t)εR^(m) is the input vector, y(t)εR^(p) is the output vector, and M,NεR^(nxn), BεR^(nxm), and DεR^(nxp) are so-called the MNA matrices. M and N containing capacitances in C, inductances in L, conductances in G and resistances in R are positive definite; B. let A=−(N+s₀M)⁻¹M and R=(N+s₀M)⁻¹B, where N+s₀M is assumed nonsingular, the Taylor series expansion of X(s) about so being given by ${{X(s)} = {\sum\limits_{i = 0}^{\infty}{{X^{(i)}\left( s_{0} \right)}\left( {s - s_{0}} \right)}}},{where}$ X^((i))(s₀) = A^(i)R is the ith-order system moment about s₀; and C. the notation V_(q(L)) and W_(q(L)) are used to represent the biorthogonal bases yielded from the block Lanczos algorithm from the Krylov subspaces K(A,R,k) and K(A^(T),D,k).
 2. The generalizations of adjoint networks techniques for RLC interconnects model-order reductions as claimed in claim 1, the purpose of the moment matching being to construct a reduced-order system such that Ŷ^((i))(s₀)=Y^((i))(s₀), wherein A. the ith-order output moment about so is calculated as Y^((i))(s₀)=D^(T)X^((i))(s₀); and B. the model-order reduction problem is to seek a q-order system, where q<<n, such that ${{\hat{M}\frac{\mathbb{d}{\hat{x}(t)}}{\mathbb{d}t}} = {{{- \hat{N}}{\hat{x}(t)}} + {\hat{B}{u(t)}}}},{{\hat{y}(t)} = {{\hat{D}}^{T}{\hat{x}(t)}}}$ where {circumflex over (x)}(t)εR^(q),{circumflex over (M)},{circumflex over (N)}εR^(qxm), and {circumflex over (D)}εR^(qxp), and the corresponding ith-order system moment and output moment about s₀ is {circumflex over (X)}^((i))(s₀)=(−({circumflex over (N)}+s₀{circumflex over (M)})⁻¹{circumflex over (M)})^(i)({circumflex over (N)}+s₀{circumflex over (M)})⁻¹{circumflex over (B)} and Ŷ^((i))(s₀)={circumflex over (D)}^(T){circumflex over (X)}^((i))(s₀). 